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Nonperturbative pair production in interpolating fields 
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We compare the effects of timelike, lightlike and spacelike one-dimensional Inhomogeneities on 
the probability of nonperturbative pair production in strong fields. Using interpolating coordinates 
we give a unifying picture in which the effect of the inhomogeneity is encoded in branch cuts and 
poles circulated by complex worldline instantons. For spacelike inhomogeneities the length of the 
cut is related to the existence of critical points, while for lightlike inhomogeneities the cut contracts 
to a pole and the instantons become contractable to points, leading to simplifications particular to 
the lightlike case. We calculate the effective action in fields with up to three nonzero components, 
and investigate its behaviour under changes in field dependence. 
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I. INTRODUCTION 

Pair production in strong fields gives us insight into 
nonperturbative fundamental physics m- Advances in 
technology have in recent years driven interest in ob¬ 
serving nonperturbative pair production using intense 
lasers mn], which presents a veritable experimental 
challenge even with optimally focussed laser light m- 

It has been found that the pair production probability 
tends to be higher (relative to the locally constant ap¬ 
proximation) in time-dependent electric fields E(t), and 
lower in position-dependent inhomogeneous electric fields 
E{z) [13]. In the latter case there is even a critical point 
beyond which the probability is identically zero [niini. 
Between these two cases is an electric field with light¬ 
like inhomogeneities, E{t + z). In this case the proba¬ 
bility is given exactly by the locally constant approxima¬ 
tion nnunj. 

Here we would like to understand more about how the 
spacetime dependence of an electric field affects the pair 
production probability. Note that the three cases above 
cannot be related by a Lorentz transformation. In order 
to investigate the transition between spacelike and time¬ 
like field inhomogeneities we will therefore consider elec¬ 
tric fields depending on a single interpolating coordinate 
of the form (1 — a)t + az for a S [0,1], following 

We will use the worldline formalism (see |21l 1^ . and 
e.g. |23It2^ for recent applications), and in particular 
worldline instantons m m [iniii7|. These are peri¬ 
odic, in general complex |28fl5T| . solutions to the clas¬ 
sical equations of motion. The classical action evalu¬ 
ated on these solutions, together with the contributions 
of fluctuations around the instantons, gives a semiclassi- 
cal approximation of the effective action, and of the pair 
production probability. We will see that this contribu¬ 
tion is an integral over the instanton itself. Since the 
instantons are complex, these contributions depend, by 
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Cauchy’s integral theorem, only on the structures which 
the instantons circulate in the complex plane. 

This allows us to unify, and extend, previous obser¬ 
vations on the impact of field inhomogeneities on pair 
production. We will see that for timelike and spacelike 
inhomogeneities the instantons circulate differently orien¬ 
tated branch cuts and so are fundamentally extended ob¬ 
jects. As the field dependence becomes lightlike, though, 
the branch points coalesce and become poles, so that the 
instantons are contractable to points |32|. We will see 
that this leads to the known localisation of the effective 
action in the lightlike limit [ini[i7|. 

Another motivation for this paper is recent interest in 
pair creation in two-component, rotating fields, which 
model the electric antinodes of a circularly polarised 
standing wave |33HSS]- We will here extend these calcu¬ 
lations to electric fields depending on our interpolating 
coordinates, and with up to three components. 

This paper is organised as follows. In the remainder 
of this introductory section we recall the worldline ap¬ 
proach to pair production in strong fields, introduce our 
interpolating fields and outline the calculation to be per¬ 
formed. In Sect. |T^ we present an explicit example of the 
instantons in an interpolating Sauter pulse, and relate 
the field inhomogeneity to the structures circulated by 
the instanton. In Sect. |III| we complete the calculation 
of the effective action in electric fields with up to three 
nonzero components depending on the interpolating co¬ 
ordinate. In Sect. |IV| we analyse the behaviour of the 
effective action as a function of the interpolating param¬ 
eter in several examples. We conclude in Sect. |V| 


A. Pair creation in interpolating fields 

The probability P of pair production in an external 
field is related to the effective action T by 

Ppairs = 1 - , ( 1 ) 

Our starting point is the one-loop worldline representa¬ 
tion of r in terms of an integral over proper time T and 
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a path integral over periodic paths , 

OO 

r = /f 

0 

where the classical action is 


S = 


m^T 

2 


+ 



0 




eA^{x)x^ 


( 2 ) 


( 3 ) 


Here r parameterises the worldline. Our first aim is to 
examine how F depends on fields of given strength, shape 
and direction as the field dependence, or inhomogeneity, 
varies from temporal to spatial. To do so we introduce 
interpolating coordinates {q,d} as [ISfEO] 



To calculate F we follow essentially the same steps as 
in m- We first find the periodic instanton solutions 
of the classical equations of motion in order to identify 
the dominant, exponential, contributions to F. Once we 
have the instantons, we evaluate their classical action, 
which is the saddle point contribution to F. We then in¬ 
tegrate over fluctuations around the instantons in order 
to obtain a prefactor contribution. Finally we perform 
the T-integral, also with a saddle point approximation, 
in order to obtain the final expression for F. This cal¬ 
culation will reveal new insights into the structure and 
properties of complex instantons. Regarding this, for the 
symmetric fields {E{t) = E{—t)) usually considered, the 
exponent in Q becomes real after rotating both proper 
time T and time t to Euclidean space. However, this ro¬ 
tation does not make the exponent real for general fields, 
so the instantons will in general be complex | 28 fl!n] . We 
will therefore use Minkowski coordinates throughout. 

Our results will be valid in the semiclassical regime, 
i.e. for weak fields and not too large adiabaticity 7 , 


in which 0 G [0, tt]. This is not a Lorentz transfor¬ 
mation and so we are considering Lorentz inequivalent 
cases. We consider electric fields of a given profile E^{q) 
(i.e. given form and amplitude) which always point in 
the z-direction. As 9 varies we interpolate between time- 
dependent homogeneous electric fields E^{t) at 0 = 0, 
fields depending on lightfront time (t-\-z)l^/2 at 0 = 7r/2 
and finally spatially inhomogeneous static electric fields 
E^{z) at 0 = TT. (The coordinate d interpolates between 
position z, lightfront position {—t + z)fV2, and time —t.) 

Let A||(g) = E^{q), then we take as gauge potential 

, (5) 

where d.x = d. It is easily checked that only F 03 = E^{q) 
is nonzero. 

We will also consider a time-dependent rotating elec¬ 
tric field under the transformation ([^. For this we use 
the gauge potential, T= { 1 ) 2 }, 


where Eq and w are typical strength and frequency scales 
of the considered field. This can be understood by rescal¬ 
ing T T/{eEo) and a: —>■ x/w in ([^, which makes the 
whole action inversely proportional to Eq and the inte¬ 
gral term inversely proportional to 7 ^. See [131 [33 for 
more details. 

We now set m = 1 and absorb factors of e into the 
gauge potentials, reinstating the mass only in final re¬ 
sults. Throughout we use the notation 

c:=cos 0 , s:=sin 0 . ( 10 ) 


II. COMPLEX WORLDLINE INSTANTONS IN 
INTERPOLATING FIELDS 


, ( 6 ) 

The transverse potential gives electric and magnetic fields 

A^g) = cosfA'(g) , B\q) =-sin ^e"^A'j{q) . (7) 

For 6 * = 0 (|^ describes a time-dependent, rotating elec¬ 
tric field, for 0 = 7 r /2 a plane wave, and for 0 = tt a 
static, inhomogeneous magnetic field. 

Note that the amplitude of the transverse fields trans¬ 
forms as we rotate. Nevertheless we present the calcu¬ 
lation of the effective action in the two types of field 
together, i.e. we will consider three-component electric 
fields, with the ‘longitudinal’ (strictly, z) and transverse 
components described by (§ and respectively, the 
potential being simply the sum of these, 

= A^^{q)d + A^{q)x^ . ( 8 ) 


In our interpolating coordinates the equations of mo¬ 
tion are 


cq — sd = TA'^x'^ , 

cd + sij = TA'^^q , ( 11 ) 

x-^ = TA'^q , 

which are just the Lorentz force equations. The latter 
two equations are readily solved for d and x^ in terms 
of q, 


d 


x^ 


-\T{A{q)-p)\\ - sq 
T{A{q) -p)^ , 


( 12 ) 

(13) 


where the integrations constants py are determined by 
integrating (12 1 and (13); periodicity of the instantons 









3 


requires that be equal to the r-average of the poten¬ 
tial, 


Phi- = j dT^ii,-L(g) =: (^)ii,i ■ (14) 

0 


This means that p depends, in general, on c, 7 and proper 
time T |33]. We also find 


x^Xu = const =: , 


(15) 


which defines a fourth constant o, see below. Using (12 1 , 
(131 and (151 yields an expression for g as a function of q] 

{ca^ + (^ - P)y + c(^ - p)i) . ( 16 ) 


To understand what p represents it is helpful to rewrite 
the equation of motion for g(r) as 


q = T\j^,{A-p\+cA'M-p)i-) ■ ( 17 ) 


and then to compare with the phase-space equations of 
motion, see e.g. m- In our interpolating coordinates the 
phase-space equation of motion for g(r) is precisely 0 
with appearing as the canonical momenta (conju¬ 
gate to d.x and x^ respectively) of the produced electron- 
positron pair. We will show below that the values of 
required by periodicity correspond to the saddle point of 
WKB integrals over Py.^. Note that A and p appear only 
in the gauge invariant combination A — p, which obeys 


(A - p) = 0 [sa. 

The solutions to (16 1 , and to the other equations of mo¬ 
tion, are in general complex, i.e. closed curves in the com¬ 
plex plane. Taking a square root of (16) implies that the 
velocity q has a branch cut in the complex g-plane, with 
the branch points corresponding to the ‘turning points’ 
where q = 0. It is at this stage useful to consider ex¬ 
amples of the instanton solutions, before going on to the 
calculation of the effective action. 


A. The Sauter pulse 

For our example we turn off the transverse fields, and 
take for the longitudinal electric field the well-studied 
Sauter pulse 

E^{q) = Eosech^ {ujq) . (18) 



FIG. 1. The complex instantons ^(t) in a time-dependent 
sech'^ electric field, o = c = 7 = 1 and various imaginary 
tq. As Im To —>■ 0 the instantons become purely imaginary 
(purely real in Euclidean space) and sit on the branch cut, 
shown as the dashed/black line. The instantons cannot cross 
the poles in the field at ±i7r/2, shown as dots. 


where 7 is the adiabacity parameter from § , 7 = 07 
with a as in 0 , To is a constant and n G Z is the turning 
number of the instanton. The constant a is related to 
proper time T via the relation 


cEqT^/i -b 7^0 = —2mTi . 


( 21 ) 


As a check, we note that ([19[)-(21 \ agree with the solu¬ 
tions in P7] in the limiiQ c —>■ 1 . 

Observe first that a real tq can be absorbed into a 
redefinition (reparameterisation) of t. In this case we 
obtain a purely imaginary q^r) and a complex d{T). The 
instanton ^(t) oscillates back and forth along a straight 
line between the two turning /branch points where q van¬ 
ishes; using (19) we have: 


• 2 rTi2 1 2 

q = 1 I ca 


1 


tanh^ {ujq) 


( 22 ) 


The shape and strength of this field do not change under 
the transformation so here we examine purely the 
dependence on the spacetime inhomogeneity. 

The instanton solutions to (11) are 


q{T)= — sinh 


i'f^/c 


\/U 


■ 7^0 


d{T)= - q{T) -b 

c 


cw 1/1 -b 7^0 


sinh 


-1 


sin 2 n 7 r(T — tq) 


jy/ccos 2 n 7 r(r — tq) 


(19) 

( 20 ) 


Consider then an imaginary tq. This gives a com¬ 
plex instanton q^r) that forms a loop around the turning 
points in the complex q plane |32l [37H^ . As the imagi¬ 
nary part of Tq becomes larger, the size of the instanton 
loop increases until it encounters poles in the potential at 
(jjq = EiT 12. There are no periodic solutions beyond this 


^ Our notation differs from that in as both t and T are there 
rotated to Euclidean space. Also our T is a factor of 2 larger. 
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FIG. 2. The complex instantons qir) from Fig. plotted 
at a height of Im uq using (221, in order to to show how they 
circle the branch cut in the instanton velocity q. 


1.0 - Re[ajf] 


/ \ - lm[a; f] 
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point; similar behaviour is seen in dynamically assisted 
pair production schemes HiOlllI]. 

In Fig. [^we plot the instantons qir) for c = 7 = 1, 
various rg, and a = 1 ; the latter gives, as we will see be¬ 
low, the dominant contribution to F. Both the small and 
large (imaginary part of) rg behaviours described above 
can be seen. In Fig. we plot the same instantons but in 
three dimensions, using Im q as the third dimension; this 
shows how the instantons circulate the branch cut in the 
velocity, and contract around it as the imaginary part of 
Tg goes to zero. 

These complex instantons generalise the Euclidean-real 
solutions in m- Very similar structures to those found 
here are seen in complexified classical motion |42| and in 
the WKB /phase-integral formalism |87l [38] . We remark 
that in the case of time-dependent fields, i.e. c —>■ 1 , then 
q —>■ t is purely imaginary for rg = 0 , while d —>■ z is 
real. This would imply an entirely real instanton after 
rotating t to Euclidean space. For general c, though, we 
have that 


t = q cos I — d sin I , z = qsin^ + d cos ! 


(23) 


so that even for rg = 0 both t and z will in general 
have both real and imaginary parts, even after rotating 
to Euclidean space. See Fig. [^for an example. 

We now turn to the instanton behaviour as a function 
of c. A convenient way to visualise the transition be¬ 
tween timelike and spacelike inhomogeneities is to plot 
the streamlines of the instanton velocity q, using ( 221 , as 


a vector field in the complex g-plane. (This method is 
well suited for studying instantons in more general field 
shapes as we do not need to solve the equations of motion 


FIG. 3. The instantons t(r) and z(r) for c = 1/2, i.e. an 
electric field E{t + z\/3), and y — a = 1. Even for rg = 0 the 
instantons have both real and Imaginary parts. 


to obtain the stream plots.) Fig.|^shows the streamlines 
for C 7 ^ = 3, 0.1,0, —0.1, —0.5. The branch cut which con¬ 
nects the turning points <7 = 0 and which is encircled by 
the worldline instantons is highlighted. For each value of 
C 7 ^, the largest possible instanton is that which brushes 
the poles in the potential; the streamlines to the left and 
right do not form periodic loopf|^ 

For timelike inhomogeneities, c > 0, the branch cut 
lies along the imaginary axis. From (19l-(20l, periodic 
instanton solutions exist for 7 arbitrarily large or small. 
For larger 07 ^ the turning points approach the poles, but 
recede from them as 07 ^ decreases, and the length of the 
branch cut decreases. 

For c < 0 the field inhomogeneity is spacelike, and 
the branch cut circulated by q{T) lies along the real axis. 
The cut grows without limit as — 07 ^ increases, and at 
—C 7 ^ = 1 extends all the way to infinity and cannot be 
circulated. There can be no per iodic solutions for — 07 ^ > 
1 as can be confirmed from (191: we must have \/|c [7 < 1 
as otherwise the argument mThe square brackets in (191 
hits [independent of the value of rg) the branch cut of 
sinh“^ and the instanton fails to be periodic. 


^ As the chosen field is periodic in the imaginary direction, so too 
are the streamlines; parts of the instantons in adjacent periods 
can be seen at the upper and lower edges of the plots, separated 
by the poles. 
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FIG. 4. Streamplots in the complex q-plane for the Sauter 
pulse instantons. Top to bottom: C 7 ^ = 3,0.1, 0,—0.1,—0.5. 
Arrows show the direction of the velocity q. 



FIG. 5. Complex instantons for c = —1/2, i.e. an electric 
field depending on a spacelike coordinate z). A fixed 

To = i/300 is chosen, and 7 varies up to the cutoff value 
of 7 = \/2 from (241, beyond which there are no periodic 
solutions. 


This parameter constraint is, again for a = 1 and rein¬ 
serting the electron mass, 



(24) 


Hence an electric field with spatial inhomogeneity and 
given strength must have a certain frequency, or width, 
in order to be able to produce pairs. This ‘criticality’ 
is already well-known and has been studied explicitly for 
the case c = — 1 mini [27]. Here we see the simple gener¬ 
alisation to the case that the field depends on a spacelike 
combination of t and z. Note that the manifestation of 
this physics in the worldline formalism is not that the in¬ 
stantons fail to be real for > 1 , but that they fail 

to be periodic. This is demonstrated explicitly in Fig. 
in which we plot the instantons for the field E(t^/3 + z). 

We will find below that the constraint must also be 
satisfied in order to have nonzero Im T - for discussions 
of the connection between the existence of instantons and 
the possibility of pair production, see EHilSl. 

As c increases from negative values toward 0, i.e. as 
the field inhomogeneities become ‘less’ spacelike, the con¬ 
straint on 7 weakens. The constraint vanishes precisely 
when c = 0 , the point at which the inhomogeneities be¬ 
come lightlike, and beyond which they are timelike again. 

Fig.jH shows that as c —>■ 0 , from either above or be¬ 
low, the turning points coalesce and the branch cut dis¬ 
appears. What remains is a (simple) pole; the instan¬ 
tons in the lightlike case are therefore contractable to 
points, whereas the instantons in all other cases c yf 0 
are contractable only around extended objects, namely 
the branch cuts. This is what singles out the case of light¬ 
like inhomogeneities. It is in just this case, where we lose 
the extended structure within the instanton loops, that 
the effective action T agrees with the locally constant 
field approximation [ISIlIIllll. That the instantons are 
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contractable to points is consistent with the direct calcu¬ 
lation of ([^ in |35], in which the Minkowski space loops 
contributing to T were seen to localise in surfaces of con¬ 
stant lightfront time. 

To recover a given instanton solution of the lightlike 
case c = 0 from the instantons with c ^ 0 we must scale 
some of the instanton parameters with c. Consider (19) 
in the limit c —>■ 0. If we want to obtain a finite sized 
q in this limit then we must scale Tq such that either 


{e 


,i2-n-To 


00 , e 


— i27rTo 


0} or vice versa. The limit of 


(191 is then 


9 =-sinh-i , (25) 

CJ 

in which r is a constant; these are indeed the solutions 
found in [32] • We will show below that we do not need to 
perform any extra rescaling in the final result for the 
pair production probability, i.e. the observable, which 
will have a smooth limit as c —> 0. 


III. CALCULATION OF THE EFFECTIVE 
ACTION 


We proceed to calculate T in the combined fields ([^ 
and describing three-component electric fields and 
two-component magnetic fields depending on q. We allow 
for more general field shapes than the symmetric fields 
(E(q) = E(—q)) previously considered. 


A. Classical action 

Using the equations of motion, the classical action of 
an instanton solution is 

1 1 

5'= + ^ J dr ^ JqiA^^-p^^) . (26) 

0 0 

Now we observe that all integrals over r G [0,1] can be 
expressed as complex contour integrals over the instan¬ 
tons themselves by trading dr for dq. Consider the final 
term in the above which vanishes for the commonly stud¬ 
ied cases of t-dependent or z-dependent fields, i.e. c = ±1 
and s = 0. If q were real we would write this term as a to¬ 
tal derivative and drop it due to periodicity. If we instead 
write it as a contour integral over g(r) it becomes 


introduce the function G to compactify notation: 


G(a?,p) 


J dTqyicT^) 

0 


y* dq^Jca‘^ {A — p)2 + c{A - p) 


2 

_L ’ 


(28) 


where the same branch is chosen as for the velocity, 
from (161. The remainder of the action can then be writ¬ 
ten 


S=^{l-a^) + G{a^,p) . (29) 

As the classical action is an integral over the complex 
instanton itself, it must be invariant under contour de¬ 
formation of the instanton, as the integrand in G is an¬ 
alytic away from the branch cut and poles. Changes in 
the parameter tq introduced in Sect.[II|simply correspond 
to those deformations for which the instanton remains a 
solution to the classical equations of motion (32). Hence 
instantons with different tq give the same contribution to 
the classical action; this is a complex generalisation of the 
reparameterisation invariance which allows us to choose 
an arbitrary real tq. This is useful, see below, in defining 
integrals which would otherwise apparently be divergent. 
This extends the results of [32] to both time-dependent 
and position dependent electric fields. 

We will find that all terms in the effective action T can 
be expressed in terms of G and and its derivatives 


Go := da^G Gi := dp^G etc., (30) 

where i = {||, T}. G is closely related to g in except 
that we include more arguments corresponding to the 
momenta p. Essentially the same function appears in the 
WKB treatment of pair production in |85| . 

Changing variables from r to q, we can also write down 
two periodicity constraints 

1 = (1) = ^ ^ , (31) 

Pi = (A,) ^ Q = <j> dg^^wJ^ ^ (32) 

which (implicitly) determine a and p, as previously found 
for real instantons in various fields |33l In 

terms of G these constraints become 


l = |Go(a^(r),p(T)), 

(33) 

0 = G,{a\T),p{T)) . 

(34) 


- ^^dg (A||(g)-p||) , (27) 

which will indeed vanish provided the instanton does not 
circulate poles of the potential. We will assume this in 
what follows. Turning to the remainder of the action, we 


B. The fluctuation determinant 

We now turn to the contribution of fluctuations about 
the instantons found above. We will be brief and high¬ 
light only differences between this and existing, similar, 
calculations in the literature. 
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The second variation of the action in ([^ with respect 
to a fluctuation 5x around an instanton solution is 


1 

<S '2 = y ^ {5q Sd 6x^ Sx'^) A 


where 



(35) 


It is now straightforward to impose the boundary con¬ 
ditions (381, which leads to long and not particularly 
illuminating expressions for For r = 1 though all 

these expressions can be written in terms of derivatives 
of G. We encounter for example terms such as 


1 = -—G 

cT3 ' 


00 


(42) 


A = 


sd'^ 

0 0 1 



0 

0 ^2 0 



V 0 

0 0 52/ 



IH 

c<i-|-sg\ 1 f)( ] 

q ) ^ q \ q ) 

cd+sq g g 

q q ^ 

q 


Q cd+sq 

6 

0 0 

0 


-d^ 

d 

0 0 

0 

V 

-dt 

0 0 

0 / 


(36) 

The Sx integral is Gaussian and gives the determinant 
of A, which can be computed with the Gelfand-Yaglom 
method as in m- Zero modes are avoided by means of 
Dirichlet boundary conditions, Sx^(O) = Sx^(l) = 0, and 
an integral over the initial instanton position x^(0). 

The determinant is given bj0 detA = det“i^^(l), 
where °‘(f> are the four solutions of the Jacobi equation 

A(/> = 0 (37) 

satisfying 

“,^(0)=0 = a,/3 = 1,2,3,4 , (38) 

where the upper left indices distinguish the four solutions 
while the upper right denote yect or components. We And 
the solutions by multiplying (371 with [q d x^ x^) and 
(O Sij ), j = 1, 2, 3, and expanding in terms of the trivial 
solutions, 


= °‘h 




( “ \ 

d 

+ 

add 

x^ 


\P) 




(39) 


Solutions can be expressed in terms of the constants 

«^l,2,3.4 

= 4 ( - -f °‘k^{sq + ch) + “fc-"cx^) , (40) 


and 


= -{°‘e -°‘h{sq + cd)) 


^d^ = - °‘hx^ 


where the left hand side appears, at first sight, divergent 
because of the turning points <7 = 0. However, these sin¬ 
gularities are avoided by the complex instantons which 
circulate the turning points, so that q ^ 0 along the 
whole contour. We can then, if we choose, deform the 
contour in G down to the branch cut without encounter¬ 
ing any divergences. Turning point singularities are also 
avoided with such contours in the WKB/phase-integral 
formalism I371E9]. After some algebra we And that the 
determinant can be compactly written as|^ 


det A = det“())^(l) 

= -^(^00 - Go^G-'Go,) det G.,- 


(43) 


where i = {|j, T}. This reduces to the determinant in |35| 
for fields with Goi = 0 (e.g. for p = 0) and c = 1, and 
to the determinant in m for one-component fields with 
c = ±1 and p = 0. In general though Goi is nonzero. 
This is in a sense a moot point, though, as everything 
in the round brackets in (431 will soon cancel against a 


similar term coming from the proper time integral. 

We turn to the x'^(O) integrals; those over d(0) and 
Xi(0) give volume factors. Including the square root of 
the determinant in (431 we have the integral 




(44) 


If the go-contour were the same as in (331, then (441 
would equal one. Instead (441 contributes a factor of lj2. 
This is shown for real Euclidean instantons in ^7}. To 
see it in the complex case, we propose a contour that 
starts (ends) on the real axis to the left (right) of the 
branch cut, and passes either above or below it. In the 
lightfront limit, the branch cut shrinks to a pole and the 
usual ie-prescription, rn? —> w? —ie in (|^, would instruct 
us whether to go above or below. As the instantons are 
symmetric under reflection about the real axis, the inte¬ 
gral (441 contributes a factor of 1/2. 

Finally we perform the proper time integral in a sad¬ 
dle point approximation. When varying the action (291 


(41) 


in order to obtain the saddle point equations, one must 


^ This method applies to the ratio of two determinants. We have 
divided by the free determinant, which is equal to one in our 
conventions. 


In our approximation we consider only the instantons with lowest 
turning number, n = 1. As a result the Morse index is fixed by 
the requirement that the probability be positive. 























remember that a and p can depend on T. Here the two 


conditions (331 and (341 for a and p become nseful and 


greatly simplify the proper time derivative: 


dT dT dT 

lil-a 


dT 




(45) 


Hence the saddle point is determined simply by setting 
= 1, as snggested earlier. At the saddle point we have 

S ^ G{l,p) , T —>• 2G'o(l,_p) , Gi{l,p)=0, (46) 

where the final equation determines p in terms of c and 
other field parameters. We also need the second variation 
of the action with respect to T. After some simplification, 


again using (331 and (341, this becomes 
d2,5 


dT 


)(< 


't-i 


= -^(Goo-Go.G-^Go, 


(47) 


in which we recognise the same factor as appears in (431, 
with = 1. 


C. Final result 

Collecting terms from the proper time and path inte¬ 
grals we find 


Im r = C 




(48) 


where 


Stt^ ^det{iGij) ’ 

G{l,p) = dg^ c+{A-p)f^+ c{A - p)l (49) 


and are determined by 


Gfyl,p)=0 


(50) 


The large round brackets in (431 and (471 have canceled, 


a simplification that might have been expected given 
Gutzwiller’s trace formula |33|- That iG is real follows 
from the fact that the contour in G can be chosen along 
an instanton, and that the integrand is the velocity up to 
a factor of i. As noted in I27| we do not need to find the 
instantons to evaluate ( |48[ ); q is now only an integration 
variable and we can choose any contour that circulates 
the branch between the turning points where the square 
root in G vanishes. 


The final result (481 may also be obtained by using the 


saddle point method for a WKB momentum integral 


d^p 


Im T = — , 

2 J (27r)3 


iG 


(51) 


It was shown in [IS] that for longitudinal fields with 
antisymmetric monotonic potentials (with c = — 1) the 


wordline instanton results in m agree with WKB /phase- 
integral results with saddle point at zero canonical mo¬ 
mentum (p = 0). Equivalence between the worldline 
formalism and WKB was also shown in [3^ for three- 
components fields with p = 0 and c = 1. Now we have 
shown that the two formalisms are equivalent also for 
more general field shapes with p fy 0 and for — 1 < c < 1. 

Consider now the transition between the cases of time- 
dependent and position-dependent fields, as c changes 
sign. In the limit c —>■ 0 the transverse fields become 
a plane wave and drop out. Further, all square roots 
drop out and the branch points and cut coalesce into a 
pole that can be used to perform the contour integrals 
in G and its derivatives using Cauchy’s residue theorem, 
as in [32]. The instantons are deformable to the poles, 
i.e. to points in this limit. It is the contraction of the 
branch cut which leads to the localisation of the instan¬ 


tons in the lightfront limit. The residue for (341 becomes 
proportional to the derivative of the electric field, imply¬ 
ing [32] 


7'ii=^ii(9) and {q) = 0 , 


(52) 


so that the instanton circles an extrema of the electric 
field. 

In the prefactor calculation, the factor in the large 
round brackets in the fluctuation determinant (43 \ van¬ 


ishes, which signals the presence of zero modes and sug¬ 
gests that, although the final results are simple on the 
lightfront, the calculation can be more subtle if it is per¬ 
formed at c = 0 from the outset. In this sense nonzero c 
acts as a regulator (which was one of the original moti¬ 
vations for introducing coordinates interpolating between 
instant-form time t and front-form t + z in |18H2()| 1 and 
allows us to complete the calculation in |32| by extend¬ 
ing the localisation seen in the instantons to the whole 
effective action. In the limit c —>■ 0 the effective action 
becomes 


Im r 




167r3 




3 t^2 . I /E 


-E" 


(53) 


where, from (52), the field is evaluated at its maximum. 


(531 agrees with the locally constant field approximation. 


IV. TIMELIKE VS. SPACELIKE 
INHOMOGENEITIES 


We now study examples of the effective action (481, 
beginning with purely longitudinal fields, = 0. 


A. Longitudinal fields 


For this case we have Gi, oc fy, and (48) reduces to 


Im T = 




r^3 


327r2 iGo^-iGm/c ' 


(54) 
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To generalise to spinor QED, we simply have to multiply 
(541 with a factor of 2, as can be shown in the same way 
as in [131121] • 

To make the dependence on c more obvious we will take 
the usual case of antisymmetric, monotonically increasing 
potentials A{q) with p = 0; in short the same field shapes 
as in m but now depending on our general coordinate 
q. Let A[q) = fiujq)/^, in which / is a dimensionless, 
monotonically increasing shape function with — 1 < / < 
1. Changing variable from q to y defined by f{ujq) = 
we find 


G(1,0) = -*^5(c7^), 

eEo 


where g is m 


5(c7^) := 



(55) 


(56) 


with f expressed in terms of p = —c^'^y^. We note 
briefly for this class of fields that, changing variables from 
g to ^ as in [nisi], one can show that 


27ri 

W)' 


(57) 


which is the same relation as was found in |82j for 
lightfront-time dependent fields (c = 0). In the case 
of constant fields, for which (eE) —?> gEq, (57) gives 
the location of the poles in the proper-time integral, the 
residues from which generate the imaginary part of the 
effective action, i.e. give pair production |49j . 

Returning to the effective action, we substitute (56) 
into (ImI to obtain 


Im r 


V3V2{eEof/‘^ 

32^" 7 5c[C5(C)]y-5|[Cff(C)] 


C=C7" 


(58) 

This generalises the result in m-- the difference is that 
a factor of c now multiplies the adiabaticity parameter 
squared, 7 ^. In m it was noted that the result for A{x) 
is obtained from that of A{t) by 7 ^ —)■ — 7 ^, which in our 
interpolating coordinates corresponds to taking c from 1 
all the way to —1. Importantly, the effective action is a 
smooth function for all c, and in particular is continuous 
as c goes through 0 . 

To illustrate, take the Sauter pulse ( |I^ . Reinstating 
the electron mass and writing E = eEQ/rrr = Eq/E g, the 
ratio of the peak field to the Schwinger field, one finds 


Im r 


£3/2 

IbTT^y 


(1 + c7^)3/4 


exp 


TT 2 \ 

1 -I- a/TTc^/ 


(59) 

see also [27] and, for higher-order corrections, |3^Fig.[^ 
shows a contour plot of the effective action ( [5^ as a 
function of c and fieldstrength £ at fixed frequency. We 



FIG. 6. Contour plot of the effective action Im F /V^ as in 
(591, in units of the electron mass. In order to show struc¬ 
ture, the frequency is fixed to the high value of w = m/2. For 
c < 0, when the field inhomogeneity is spacelike, the electric 
field strength must obey Eo/Es > \/|cf/2 in order to be ca¬ 
pable of producing pairs. We see that to achieve a given pair 
production probability it is preferable for the field to have 
temporal rather than spatial inhomogeneities. 


see that the probability increases from zero beyond the 
critical curve 


C 7 ^ = —1 . 


(60) 


This curve connects the known critical point 7 = 1 for 
the case of spatially inhomogeneous fields E{z) to j = 00 
at c = 0 , where the field dependency becomes lightlike, 
E{z) — E{t + z). For c > 0 on the other hand, in 
which the field dependency is time like, there are no criti¬ 
cal points. The same behaviour is seen in Fig. where we 
plot the effective action as a function of c and frequency 
w/m at fixed field strength. Hence the lightfront limit 
c = 0 separates systems with critical points from those 
without. It is already known that for q — t the pair 
production probability is larger than given by the locally 
constant approximation, while for g = 2 ; it is smaller, and 
for the lightlike case, q = x*, equal; the contours in Fig.[^ 
confirm that (at fixed field strength and frequency) any 
timelike dependence yields a higher pair production rate 
than any spacelike dependence. 

It is worth saying that for the near future’s physically 
realisable parameters, i.e. electric fields delivered by in¬ 
tense laser systems of optical frequency u! ^ leV and field 
strengths of order at most Eg/lOO, (59) varies only im¬ 
perceptibly, at the order of 1 part in IO 7 as we interpolate 
between spacelike and timelike inhomogeneities. Hence 
the biggest distinction between the spacelike and timelike 
cases comes from the different volume factors appearing 
which, after regularisation, go roughly like the Rayleigh 
range for c = -1-1 and the pulse length for c = — 1. (The 
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FIG. 7. Contour plot of the effective action Im F/F® as in 
( |59| l, in nnits of the electron mass. Here E = Es/W, and the 
critical curve is at 100c(a;/m)^ = — 1. 


FIG. 8. Effect of a transverse field on the criticality condi¬ 
tion, compare ( |64[ | with ( |60[ ). Parameters as in Fig. For 
fields with spacelike inhomogeneities c < 0, turning on trans¬ 
verse electric and magnetic fields reduces the parameter region 
in which pair production is possible. 


remaining area factor, common to both volumes, is the 
focal spot size.) 

There is a simple argument which allows us to 
check (58l. Let us make the frequency/wave length ex¬ 
plicit in the argument of the field, writing E = E(ujq). 
For, e.g. 0 < 6 < 7 r /2 we boost with velocity v = 
tan( 0 / 2 ), parallel to the field, such that the interpolat¬ 
ing coordinate q becomes proportional to time in the new 
frame, q = \fct'. In this new frame we have a time depen¬ 
dent electric field E[^s/ajjt') HTHSHIEI], with frequency 
scaling with c as 


uJ = \fcu: . 


(61) 


Further, the volume in the d-direction is (for |f'| 

V) related to the volume in the x'-direction by V' = 
yJcV^. Hence, by partially undoing our rotation with a 
Lorentz transformation, we can check (541 by rescaling 
parameters in the result of m- 


For c < 0 we can similarly boost to a frame with a z' 
dependent field and frequency (or rather wave vector) 
scaling as w' = y/^oj. If there are transverse compo¬ 
nents the same Lorentz transformation simplifies them: 
for c > 0 (c < 0 ) the magnetic (electric) components van¬ 
ish after the boost. However, the magnitude of the trans¬ 
verse components also scales with i/±c, which means 
that the corresponding adiabaticity parameters are in¬ 
dependent of c. 


This argument also informs the lightlike limit c —0. 
In the final expression for T we can see that c —>■ 0 takes 
the adiabacity to zero and, comparing with (611, is ef¬ 
fectively a zero frequency limit. This is again consistent 
with the lightlike case agreeing with the locally constant 
field approximation. 


B. Longitudinal vs. transverse components 


Consider now adding an additional electric field com¬ 
ponent which points in the transverse direction (along 
with the magnetic component) as in ([^. We take Ai 7 ^ 0. 
Let this and the longitudinal field component have the 
same shape but different magnitudes, 


^ 1 ( 9 ) —'^^ 11 ( 9 ) j (62) 

with A dimensionless. Assume again that Ay = //y is 
antisymmetric and monotonic, then we have 

G = - j>dq c+ cA2)A2 , (63) 


which, comparing with the purely longitudinal case, sim¬ 
ply amounts to rescaling 7 ^ by factor 1/(1 -|- cA^). The 
criticality condition then becomes 


cy 


1-f cA2 


> -1, 


(64) 


provided 1 -|- cA^ > 0 is also obeyed (as otherwise the 
invariant E^ —becomes negative and, E.B=0, there is 
no pair production.) The constraint (64) reduces to (60) 
for purely longitudinal fields, A = 0; the two constraints 
are compared in Fig. For c < 0 we see that the effect of 
the transverse fields is to reduce the volume of parameter 
space in which pair production is possible. 


C. Two-component ‘rotating’ fields 

We turn now to a two-component transverse field 

^^( 9 ) = —{coswg,sinwg} , (65) 
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FIG. 9. Complex instantons ojt, ujxi and ujX 2 for c = 7 = 1, 
with iot = 1 and initial condition iut{0) = —1.75. 


which for c = 1 is a time-dependent rotating field, and 
for c = 0 a transverse (monochromatic) plane wave. (For 
c < 0 one can boost to a frame where the field becomes 
purely magnetic and there is no pair production.) For 
the time-dependent case the effective action was studied 
in [Ml SSI 1^ using WKB, in [M] using the Wigner for¬ 
malism, and the Euclidean instantons found in |33| . We 
will therefore be brief here, only using this as an example 
to show that the complex structures above exist in other 
fields. We believe though that this is the first calcula¬ 
tion of F including the prefactor contribution, using the 
worldline approach. 

We begin with the instantons. Fulfilling the periodicity 
constraint (321 or (50) is essential here. Given the form 


of the field (65), we make the following ansatz for p^: 


P 


p± = —{coswg,sinwg} , 
7 


( 66 ) 


where p := 7 |pj_| is the magnitude of the transverse mo¬ 
mentum, and q is some average which, as we will shortly 
confirm, may be freely chosen. The constraint which de¬ 
termines p as a function of 7 ^ is given in |33| . This con¬ 
straint may be rewritten in terms of the complete elliptic 
integrals K and E (see |521 §17.3]) as 


0 = 




(l + p)K(-a2)_2E(-a2) 


(67) 


in which tr^ := ( 7 ^ + {p — l)^)/(4p). Numerical solution 
of (671, or estimation using the approximations in |33| , 


shows that p > 1 , which is larger than would be possible 
if the instanton were real. To see the reason for this let 
q = 0 and deform the instanton so that q{T) follows a 
straight line between the turning points, recall (16). The 


line is along the imaginary axis, so while one component 
of the average is zero {sin ujq) = 0 the other is {cos ujq) = 
(coshiwg) > 1 . 

In |33| further conditions were imposed on the instan¬ 
tons to ensure reality, but these are not needed. This is 
confirmed in Fig. where the instantons for c = 1 are 
found by numerical integration of (11). A real t and a 


real staring point t( 0 ) are chosen, and indeed the solu¬ 
tions are still complex. We have confirmed that periodic 


solutions are only found if (67) is fulfilled. While ^(t) 


remains a simple closed curve, the instantons x^{t) can 
self-intersect. 

The classical action of an instanton, G, can be written 


where h again depends on elliptic function^ 


( 68 ) 


^(7^P) = ->/7^ + {P+^Y 


K-E 


7^ + (p-l)^ 
7^ + (P + 1)^ 


)> 


(69) 

as also appears in the WKB calculation of |M|- Due 
to the symmetry of the field G does not depend on q. 
In the calculation of the effective action this leads to a 


zero mode, and the determinant of Gij appearing in (48) 


vanishes. To separate out the zero mode into a volume 
factor we follow |531 §39.4]. First write the determinant 
as a momentum integral. 




1 


det iGij ^—%G, 


33 




as holds to lowest order in the semiclassical regime, and 
change variables from p^ to p and wq; the p integral 
is performed with the saddle point method and the q 
integral gives a volume factor ixVq, or 27r per period of 
the rotating field. The result is 




-y/— det iGij 


i VtE p 

7 yJiGzsiGpp 


and the effective action becomes 


Im F 


V4c{eEoYp 
%'K^y/2hih22 ^ 


TT h\ 
ey/cEol'^J ' 


(71) 


(72) 


extending the results of imiiQii to c < 1. The dependence 
of F on c is plotted in Fig. [T^ 


V. CONCLUSIONS 

We have used a coordinate rotation to investigate 
Schwinger production in backgrounds which interpolate 
between time-dependent, homogeneous electric fields and 
inhomogeneous, static electric fields. This allowed us 
to examine the transition between Lorentz-inequivalent 
spacetime dependencies. For all field dependencies we 
found that the instanton contribution to the effective ac¬ 
tion was given by a complex contour integral over the in¬ 
stanton itself, with the physics of pair production being 
encoded in the branch cuts circulated by the instantons. 


® It is interesting to compare l |69| with Eq. 3.51 in EZ], which 
gives 5 ( 7 ^) for a longitudinal oscillating field, Eq cosa;t; the only 
difference, after using various elliptic function identities, is the 
factors of p, which are absent in the longitudinal case, and an 
overall factor of 2. 
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FIG. 10. Contour plot of the effective action (721 with 
Eo = Es- There is no pair production for c < 0. 


Note that the existence of critical points, beyond which 
there is no pair production, is not related to the reality of 
the instantons. Instead we have seen that critical points 
arise simply when the instantons fail to be periodic. 

Being complex contours, the instantons can be freely 
deformed (Cauchy’s integral theorem) around the branch 
cuts, without changing their contribution to the effective 
action. A striking property of the instantons is that they 
make this symmetry manifest: the freedom to choose tq 
in the instanton solutions represents those deformations 
for which the instanton remains a solution to the equa¬ 
tions of motion. This was previously found for the case 
of lightlike field dependencies [32] , but now we have seen 
that it holds more generally, for both spacelike and time¬ 
like dependencies, and for a range of field configurations 


including longitudinal fields and two-component rotating 
fields. 

For all timelike and spacelike inhomogeneities, the in¬ 
stantons are deformable only down to a branch cut, and 
are therefore fundamentally extended objects. We found 
though that the limit of lightlike coordinate dependence 
corresponded to a vanishing field frequency scale, in 
which the branch circulated by the instantons contracted 
to a pole. The instantons in this case are contractable to 
points, i.e. are equivalent to pointlike objects, and it is in 
just this limit that the effective action is given simply by 
the locally constant approximation. We have recovered 
this result by rotating our field dependence from time¬ 
like to spacelike, and the effective action is a continuous 
function of the interpolating parameter. 

The freedom to deform the instantons, even away from 
solutions to the equations of motion, seems consistent 
with the observation in m that the explicit form of the 
instantons is not strictly needed in order to calculate the 
semiclassical approximation to the pair production prob¬ 
ability. For the example of the Sauter pulse we have seen 
that, at fixed field strength and frequency, the probabil¬ 
ity is higher for any timelike dependence than for any 
spacelike dependence. 

We have considered fields with one dominant max¬ 
imum and consequently negligible interference effects. 
Such effects have been studied in m using a phase-space 
worldline approach, and it would be interesting to study 
the cases covered there using our formalism. 
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